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ABSTPACT 

The Plasma Spectrophotometer constructed at the Naval Postgraduate 
School in 1968 was modified by the addition of a photomultiplier cooling 
device to improve time resolution and increase the accuracy of the 
computed values of the relative intensities of the nitrogen spectrum 
in the magnetically confined nitrogen plasma column. 

The relative intensities determined were used to analyze the vari- 
ation in local electron temperature of the singly and doubly ionized 
nitrogen atoms with respect to magnetic field. The local electron temp- 
erature was found to vary with increasing magnetic field, reaching a peak 
temperature of 8.1 electron volts at 4200 gauss. 

Evaluation of the relative number densities of the two ionic 
species indicated that the local electron temperature of the column 
was the local electron temperature of the singly ionized nitrogen atom. 
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I . INTPDDUCTION 



This investigation is a continuation of the analyses of the char- 
acteristics of the highly ionized plasma column produced at the Naval 
Postgraduate School Plasma Facility which have been undertaken since 
its fabrication in 1962. The principle objective of this study was 
the analysis of shock waves in the magnetically confined plasma column 
by means of the dual -Ebert type spectrophotometer constructed at the 
Plasma Facility in 1968 by Orlicki [1] . The spectrophotometer was to 
be modified to improve time and spatial resolution and then used to 
determine the variation of electron temperatures and relative number 
densities of key ionic species with changes in the machine parameters, 
such as magnetic field. 

The key ionic species of the unperturbed plasma were identified as 
the singly and doubly ionized atom (N+ and N++) and the singly ionized 
molecule (N 2 _) [1] . These species are the basis of the continued study. 
The identified lines are listed in Table I. 



Table I. Identified Lines of the Nitrogen Plasma 



o 

N+ Wavelength (A) 



N++ Wavelength (A) 
4097.309 
4103.370 
4200.000 



o 



4124.100 

4133.654 

4145.759 

4601.490 

4607.166 



4514.900 

4518.200 



4613.887 

4621.405 

4643.106 



4535.100 
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A complete reanalysis of the plasma in the unperturbed state had 
to be accomplished first because of several modifications to the Plasma 
Facility, including the addition of a ten inch diffusion pump to the 
cathode end of the system. Thus the local electron temperature and 
the relative number densities of the two kncwn nitrogen atom species 
and their variance with respect to changes in the magnetic field were 
to be determined to identify changes brought about through incorporation 
of the modifications. 
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II. BACKGROUND 



A. SPECTROSCOPIC DIAGNOSTICS 

The nitrogen plasma produced at the Naval Postgraduate School 
Plasma Facility emits sufficient light for spectroscopic observation 
and is also sufficiently dense and long-lived so that the assumption 
can be made that the velocity of the electrons is very nearly Max- 
wellian at any instant of time and at any point in space. The local 
electron temperature can then be treated as a function of time. Since 
electron impacts dominate almost all collisional rate processes, the 
local electron temperature is essentially the quantity determined in 
spectroscopic temperature measurements. 

The electron density is usually involved in the relations used to 
calculate the electron temperature frcm observed intensities, and 
therefore temperature and density measurements are not really inde- 
pendent. An attempt to find a self-consistent solution of the various 
relationships between the directly measured spectroscopic quantities 
and densities and temperatures was made in this investigation. 

Since the assumption that the plasma is in local thermal equilibrium 
(LTE) has been made, the relationship based on the fact that densities 
in various excited states are proportional to the products of the stat- 
istical weights with the exponentials of the negative ratios of excit- 
ation energy and the thermal energy kT may be used [2] . 

Quantitatively , 




3 3 

£n(IA g'f'/I'A' gf) 



where the quantities are defined as follows: 
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1 . 



First Reference Line 



I = relative intensity 
A = wavelength in meters 

g = statistical weight of the lower state (2J + 1) 
f = absorption oscillator strength 
E = excitation potential in joules 

2. Second Reference Line — same as above except designated by prime. 

Excitation potentials and wavelengths nay be obtained frcm the multi- 
plet tables of Moore (1959) [3] , or Striganov and Sventitski (1968) [4] . 
Precise wavelengths and statistical weights can be found in almost all 
cases, but this is not true of absorption oscillator strengths. Except 
for one-electron systems, the oscillator strengths given in the litera- 
ture are approximate, with estimated errors ranging frcm 10 percent in 
the case of neutral and other relatively simple spectra to uncertainties 
of factors of two or three for more complex spectra [2,5] . 

Relative intensities can be measured with an accuracy of better than 
10 percent, unless the wavelength separation of the lines is vary large. 
Errors due to uncertainties in the oscillator strengths are usually the 
greatest source of error, and uncertainties of 30 percent in the argu- 
ment of the logarithm are not unccmmon. The uncertainty above is 
reflected in a similar temperature error, and since kT is usually of the 
order of the largest separation between excitation energies of non reso- 
nance lines, one can hardly improve accuracy by using resonance lines. 

The gain in separation between excitation energies is generally offset 
by larger errors in oscillator strengths, in relative intensity measure- 
ments, and with difficulties with self absorption [2,6]. 

Applying the concept of relative intensity to different ionic species 
to the equation for absolute intensity of a spectral line given by 
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Robinson and Lenn [6] , and using the electron temperature determined 
from the equation on page 11, one can arrive at the relative density 
equation: 



where the quantities are as defined earlier, except that they are now 
for two different species. Uncertainties from the only approximate 
validity of these equilibrium relations, frcm experimental errors, and 
also from errors in oscillator strengths, as mentioned above, are freq- 
uently of similar magnitude which implies that even with utmost care it 
is often not possible to reduce errors in relative densities below 30 
percent . 

It must be emphasized that several important criteria are applicable 
when using the above equations. These expressions are only applicable 
when the arc pressure is held constant [2] , and certain ratios of para- 
meters in both equations must be limited if meaningful results are to 
be obtained [2,6,7] . The intensity ratio must be on the order of 1.0 
to 6.0, the difference in excitation potentials must be greater than 3.0 
volts, and the wavelength separation of the lines must be small. 

It has been assumed that the electron density is sufficient to 
maintain a Maxwellian Distribution of quantum state population, and if 
this is not so, and radiative depopulation of lower energy states takes 
place to any significant degree, then the population of states having 
energies above kT would be favored and the excitation temperature indi- 
cated by the line intensities would be higher than the actual electron 
temperatures [2] . 
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B. THE NITROGEN ATOM 



The nitrogen atcm has seven identified ionization states. However, 
for this investigation, only the singly and doubly ionized states are 
of concern. Their ionization potentials are 29.593 and 47.426 electron 
volts respectively [8] . The lines which were chosen to study are listed 
in Table II along with their relative intensities and transitions. All 
of the transitions listed are allowed -[3,4] . Four lines of the singly 
ionized and three lines of the doubly ionized atoms, all marked by 
asterisk (*) , are listed among the more intense lines known in the 
nitrogen spectrum [4] . 

Table II. Relative Intensities and Transitions of 
Applicable Nitrogen Lines 



o 

Wavelength (A) 


I 


Transition 


4124.100 


4 


o 5_ 0 5 o 0 

3s P - 3p S 


4133.654 


5 


II 


4145.759 


6 


ll 


*4601.490 


11 


3s 3 P° - 3p 3 P 


*4607.166 


10 


II 


4613.887 


9 


II 


*4621.405 


10 


II 


*4643.106 


11 


II 


*4097.310 


10 


~ 2 _ 2_.0 
3s S - 3p P 


*4103.370 


9 


II 


4200.00 


6 


3s 2 P° - 3p 2 D 


*4514.900 


7 


3s 4 P° - 3?^ 


4518.200 


3 


II 


4535.100 


2 


3p 4 S - 3d 4 P° 



The state of ionization in the plasma is due to many factors. The 
primary factor is the cathode current interaction which supplies energy 
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to thermionically emit electrons which collide with the gas column. 

The higher the cathode current, the greater number of electrons emitted 
and consequently, the greater the number of collisions the greater the 
energy exchange. The electrons removed in the ionization process may 
recollide with the cathode and cause secondary emissions which are 
capable of causing further ionizations. The cathode-anode accelerating 
potential makes further energy available to the thermionically emitted 
electrons, and so it continues. Other ionizations may occur within the 
plasma due to atomic, molecular or electron interactions which would be 
functions of other machine parameters such as magnetic field. 

C. THE PLASMA FACILITY 

The experiment is carried out in the steady state plasma column, 

13 -3 

with plasma density of the order of 10 cm , provided by a hollcw 
cathode reflex arc. This plasma is confined in a ten-foot long column 
along the axis of a six coil solenoid which produces a magnetic field of 
up to 10,000 gauss and is homogeneous to within 2.5% along the axis of 
the plasma column. 

The plasma tube is a ten foot long series of four inch diameter 
pyrex sections with two access ports available on each side. The "theta" 
pinch section is located approximately thirty inches from the floating 
anode, and the access port used for the collection of spectroscopic data 
is 62 inches from the cathode (Figures 1 & 2) . 

The magnetic field is provided by six main magnet coils, a cathode 
magnet, a first anode magnet, and a mirror magnet at the floating end. 
All magnets are mounted coaxially to the tube and all are individually 
variable . 
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The nitrogen gas enters the tube at the cathode and flows through 
several centimeters of the cathode tube of tantalum where it is ionized 
by the high current discharge which runs into the interior of the tube. 
This highly ionized gas then leaves the cathode and forms the beam. 

The cathode is self heated and thus electrons are thermionically emit- 
ted from its inside walls or by field emission across a very thin sheath. 
Electrons which collide with gas molecules before returning to the walls 
are trapped and drawn out of the cathode by the axial potential gradient 
of the arc as well as by the pressure gradient of the plasma. The col- 
lisions of the electrons with the gas molecules produce positive ions 
and electrons, radiation, and possibly sane metastable excited atans. 
Nearly all of these products are trapped in the tube and return to the 
wall. This heats the wall to temperatures around 3000°C, and results 
in thermionic emission of electrons [9] . The cathode is heated locally 
to high temperatures and when cathode current is held to 60 amperes, 
the cathode life is around ten hours. The duration of the plasma column 
in a steady state was from a minimum of a few minutes to as long as four 
hours. 

Argon is used to establish the initial plasma column, and, when it 
has stabilized, the system is switched over to nitrogen and the argon is 
turned off. The transition is readily observable in the column as it 
smoothly turns from blue to pink. Data is taken when the nitrogen beam 
is established. 

D. PLASMA SPECTRDPHOTOMETER 

The plasma spectrophotometer is a dual-ebert type photoelectric 
spectrometer. It was designed to accept a beam of light through a 
single input slit, set at 300 microns, split that beam in half by means 
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Figure 3 

Internal View of Plasma Spectrophotometer 
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of a beam splitter, and pass the resultant beams to plane gratings 
which diffract the beams, after which they are focussed by spherical 
mirrors into two independent chambers [1]. (See Figure 3.) The dif- 
fracted beams exit through separate slits and are picked up by an 
Eldorado Differential Photometer Unit and recorded. The Eldorado is 
used because it contains all the apparatus necessary to perform this 
experiment, and is highly accurate, giving a high degree of spectro- 
scopic resolution. It provides the photomultiplier heads, power supply 
and recording apparatus, and is also capable of mounting both end 
window and side window photomultiplier tubes. 

The 1P21 photomultiplier tube was used for this experiment as it 

was the most suitable for the nitrogen lines under investigation. It 

o 

is a side window tube with maximum sensitivity at 4000 A and with over- 

o 

all sensitivity from 3100-5400 A. At 1024 volts supply voltage, this 

g 

tube gives a current amplification of 2 X 10 with an equivalent noise 

-12 

level of only 6 X 10 lumens. Two other tubes are available, the 

o 

6903, head-on tube, which has a maximum sensitivity at 4000 A and over- 

o 

all sensitivity from 2000-5600 A, and the 7102, also a head-on tube, 

o 

with a maximum response sensitivity at 8000 A and overall sensitivity 

o 5 

between 5600 and 9600 A. The 6903 has a current amplification of 4 X 10 

-12 

and an equivalent noise level of 6.7 X 10 lumen. It is better suited 

for work in the ultraviolet region. The 7102 has a current amplification 

5 -10 

of 4.6 X 10 and an equivalent noise level of 7.5 X 10 lumens, and 

is better suited for work in the near infra red region [10] . Using any 

or all of the aforementioned tubes gives the capability of covering the 

o 

entire spectrum of the spectrometer, 2000-9600 A. 
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Setting the required wavelengths is accomplished through the 

use of digital counters capable of locating each individual setting, 
o 

There is 1/10 A between successive settings. The gratings are mounted 
on an expanding lever-arm screw type mechanism which turns the grating 
and the digital counter at the same time. Settings can be made and 
accurately reproduced based on a prepared calibration curve (Figure 4) . 
This curve was rechecked with a mercury arc lamp, after all modifications 
were applied, and no changes were found. 
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PLASMA SPECTROPHOTOMETER 



WAVELENGTH CALIBRATION CURVE (I) 



(Number of Turns 
vs. Wavelength) 




4000 6000 (8000) (10000) (12000) (1st) 

4000 5000 6000 (2nd) 



Figure 4 
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III. EXPERIMENTAL PROCEDURES 



A. MODIFICATIONS TO SPECTROPHOTOMETER 

Since time resolution was desired, it was necessary to modify the 
differential photometer unit through the addition of a photomultiplier 
cooling system. The cooling of the photomultiplier tubes was necessary 
to reduce the detrimental effects of dark current and external noise 
[11,12] . A cooling device was constructed which reduced the dark cur- 
rent and external noise pickup rapidly and conveniently through the 
use of liquid nitrogen. No difficulty was experienced with electrical 
short circuits, but the usual problem of window fogging had to be 
resolved. 

A schematic drawing of the cooling device is seen in Figure 5. The 
side window photomultiplier tube housing from the Eldorado unit was wrap- 
ped with a coil of copper tubing ending in the exit lens housing. The 
original aluminum lens housing was replaced with a nonconducting material, 
polyvinylchloride tubing (PVC) , which had the proper internal dimensions 
and was machinable. A gas vent tube was installed in the top of the PVC 
tubing to prevent pressure buildup inside the system, and to allow water 
vapor to escape. The entire copper coil was connected to the tube housing 
with Thermon, Heat Transfer Cement [13] . Since the tube housing was 
aluminum, welding or brazing was impractical, and the heat transfer cement 
insured effective heat conduction between the metals. The entire system 
was then covered with polystyrene for further insulation. 

In operation, liquid nitrogen is placed in a reservoir connected by 
a Tee device to the input end of both coils. The liquid nitrogen fills 
the coils to the height of the liquid in the reservoir and then begins 
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evaporating. Inrnediate cooling of the tube housing takes place through 
conduction. Because the walls of the housing cool first, water vapor 
that is present condenses upon them, thus preventing harmful conden- 
sation on the tube or its wiring. The nitrogen gas flews upward through 
the coil and fills the lens housing, cooling the photomultiplier tube and 
the lens. Excess gas escapes through the vent tube. To prevent fogging, 
a 1/16 inch coil of nichrome wire, six inches long, and wrapped in a 
fiberglass sheath to prevent radiation, was wrapped around the circumfer- 
ence of each lens, and connected to a variac. Adjustment of the variac 
as needed produces sufficient heating of the lens to prevent fogging. 

Thermocouples inserted in the base of each photomultiplier tube hous- 
ing register the temperature of the resistor bank of each tube. The photo- 
multiplier tubes are cooled sufficiently for use fifteen or twenty minutes 
after filling with liquid nitrogen. About two liters of liquid nitrogen 
is needed to cool the tubes to about -80°C, and about 1/4 liter of liquid 
nitrogen should be added each 10 minutes for a thirty minute period. Data 
can be taken for about 30 minutes after this without further cooling. The 

dark current, with 1012 volts across the tube, is reduced by a factor of 
3 

about 10 within 20 minutes of applying liquid nitrogen to the reservoir. 
There is essentially no difference in the spectral sensitivity of the tube 
at low temperatures [10] . 

B. OPERATION OF THE SPECTROPHOTOMETER 

The spectrophotometer is rolled into position opposite the port to 
be used in the Plasma Facility. The unit is raised and leveled so that 
the input optics are positioned as required [1] . The paver switch of 
the Eldorado Unit is activated and one liter of liquid nitrogen is placed 
in the cooling reservoir. The variac is turned on and adjusted to heat 
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the output lenses. After fifteen minutes have elapsed, another liter 
of liquid nitrogen is added. It takes 30 minutes for the electronic 
equipment to settle out and after that interval the high voltage 
switch is activated. A lens cover is placed over the input optics 
and the meters for both photomultipliers are calibrated [14] . After 
removal of the lens cover, data may be taken. 

The first line of interest is selected by use of the digital count- 
ers and the input optics are adjusted to give maximum intensity. The 
"B" channel photomultiplier intensity reading is used as the base 
intensity in all cases, and the "A" channel intensity reading is 
adjusted to equal "B". The method known as "leap-frogging" was used 
exclusively in this investigation, i.e., starting with the same line 
on each channel, with equal intensity on both, and moving on to the 
next, registering equal intensities on both channels before moving on 
to the next [1] . Thus all intensities end up being referenced to the 
first line to come under investigation. In this investigation, it was 
always the longest wavelength under study. 

C. EXPERIMENTAL DATA AND COMPUTATIONS 

A nitrogen plasma column is established at the Plasma Facility and 
the parameters of the system, magnetic field, arc current, arc pressure, 
and anode voltage are recorded. The magnetic field desired is selected 
and held constant while all the selected nitrogen lines are investigated. 
It takes approximately 30 minutes to evaluate 14 lines. The magnetic 
field is then changed and the parameters are checked for deviation. 

This is continued until all the lines have been evaluated at every 
magnetic field in question. For this investigation, the magnetic field 
was set at 2400, 3000, 3600, 4200, and 4800 gauss. The arc current was 
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constant at 60 amps and the anode voltage varied from 64 to 75 volts 

during separate runs. The arc pressure varied during any particular 
-4 -5 

run frcm 1.6 X 10 to 7.6 X 10 torr. This will account for sane 
variance in the results. 

The intensities of each line are recorded as shewn in Figure 6, 
and the uncorrected relative intensities (URI) are computed . The URI 
of line 1 is the reading for line 1, as both channels have been balanced 
with channel "B" being the base. The URI for the remaining lines is 
computed as follows: The URI of line 2 is to the URI of line 1 as 

the average reading for line 2 is to the reading of channel "B" (base 
channel) for line 1. Any changes in readings for the same line on two 
channels (caused by moving the gratings) must be taken into account, 
and thus the use of the average reading for line 2 above. The compu- 
tations are continued substituting line 2 for line 1 and line 3 for 
line 2 proceeding through the remainder of the lines in the same manner. 
The URI must then be corrected for spectral sensitivity of the photo- 
multiplier tubes used. (See Figure 7) . The URI multiplied by the 
correction factor gives the Corrected Relative Intensity [1] . 

A computer program was written to calculate the temperature equation 
given in the "Spectroscopic Diagnostics" portion of this paper. How- 
ever, before temperatures could be computed, the values of the maximum 
and minimum intensity ratios, the maximum wavelength separation, and 
the minimum excitation potential difference had to be determined. The 
program was adjusted so that one parameter could be varied while holding 
the others constant, and was broad enough so that the parameter varied 
could assume all values in question during one run of the computer. 
Limits on the values of parameters were established so that consistant 
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PLASMA SPECTROPHOTOMETER SAMPLE DATA 



CHANNEL "A" CHANNEL "B" 

WAVE- Ionic 



LENGTH SPECIES 
APPROX, 
o 

4643 A 
4621 " 

4613 " 

4607 " 

4601 " 

4535 " 

4518 " 

4514 " 

4200 " 

4145 " 

4133 " 

4124 " 

4103 " 

4097 " 



APPROX. 


READING 


SETTING 


X 10 X 
RANGE 


221.73 


73 


220.45 


54 


219.97 


36 


219.58 


60 


219.25 


65 


215.10 


21.5 


214.19 


55 


213.96 


10.5 


196.02 


200 


193.02 


22 


192.31 


21 


191.83 


12.5 


190.67 


6.7 


190.26 


9.3 



APPROX. 


READING 


SETTING 


X 10 X 
RANGE 


207.27 


73 


206.05 


49.5 


205.52 


39 


205.14 


56 


204.77 


90 


200.73 


19.5 


199.73 


49 


199.56 


13 


181.43 


125 


178.43 


21 


177.70 


20 


177.14 


13 


176.00 


7.5 


175.64 


7.4 



I (REL) 


I (REL 


UNCORR- 


CORR- 


ECTED 


ECTED 


73 


62.78 


51.75 


44.51 


39.20 


33.71 


58.30 


50.14 


80.69 


69.39 


18.38 


16.63 


49.01 


44.35 


11.75 


10.63 


146.91 


143.24 


37.68 


37.11 


37.78 


36.23 


23.45 


23.10 


12.81 


12.62 


14.22 


14.15 



N + 

It 

ft 

It 

ft 

N++ 

tt 

n 

tt 

N + 

it 

it 

N++ 

it 



BY: J R C PURPOSE: TEMP /DENSITY CALCULATION DATA 

MAGNETIC FIELD: 3600 G AUSS ARC CURRENT: 60 A MPS 

ANODE VOLTAGE: 75 V OLTS ARC PRESSURE: 1.6 X 10~ 4 TORR 

VIEW PORT: # 2 POSITION: Centered on Beam 



RELATIVE INTENSITY CALCULATION TECHNIQUE; Lean - Frog All Referenced 

to Channel "B" 



Figure 6 
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1P21 SPECTRAL SENSITIVITY AND RESULTANT CORRECTION FACTORS 



1P21 SPECTRAL SENSITIVITY CURVE 



Relative 

Sensitivity 




NET RELATIVE INTENSITY CORRECTION FACTORS 



Wavelength Range 


Correction Factor 


Deviation 


4000-4100 ft 


0.995 


0.005 


4100-4200 " 


0.985 


0.005 


4200-4300 " 


0.975 


0.005 


4300-4400 " 


0.965 


0.005 


4400-4500 " 


0.945 


0.015 


4500-4600 " 


0.905 


0.025 


4600-4700 " 


0.860 


0.020 


4700-4800 " 


0.805 


0.025 


4800-4900 " 


0.765 


0.025 


4900-5000 " 


0.720 


0.020 



Figure 7 
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values of temperature were obtained. Three different sets of relative 
intensities determined on three separate runs of the Plasma Facility 
were used to keep the diagnosis unbiased. Limits established for the 
parameters were maximum relative intensity of 4.0, minimum relative 
intensity ratio of 0.25, maximum wavelength separation of 500 angstrons, 
and minimum excitation potential difference of 7.0 volts. The computer 
program with all corrections included and all parametric varying devices 
removed is shown on page 36. 

A computer program to calculate relative number density ratios of 
the two nitrogen species was written and modified with the inclusion 
of the above determined parameters [1] , and is shown on page 39. The 
temperature used in this computation must be taken from the temperature 
calculation above and entered as a constant. The average temperature is 
usually chosen, however any temperature within the standard deviation of 
the mean will suffice. 
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IV. RESULTS AND CONCLUSIONS 



A. RESULTS 

1. Electron Temperatures 

Entering data similar to that shown in Figure 6 into the 
electron temperature program led to quantitative results for the 
electron temperature . The local electron temperature was found to 
vary with increasing magnetic field and is shewn in Figure 8 for the 
case of the singly ionized nitrogen a team. 

The temperature varied from 3.85 electron volts at 2400 gauss 
to a peak temperature of 8.1 electron volts at 4200 gauss. It drop- 
ped rapidly to 4.4 electron volts at 4800 gauss. The compression 
of the plasma column into a smaller volume as the magnetic field is 
increased can be physically observed from the side of the Plasma 
Facility. The density of the singly ionized atoms must then increase 
and consequently atomic interactions increase giving rise to higher 
temperatures. Above 4200 gauss, the electron temperature decreases 
rapidly plausibly due to the collision rate increasing at such a rate 
as to cause a corresponding decrease in diffusion rate, and thus a 
greater distribution of energy throughout the column. These temp- 
eratures, Figure 8, closely approximate those found by Major Thomas J. 
Hay craft, USA, whose studies of the plasma column under the same con- 
ditions were conducted using probe techniques and frequency analyses. 

The temperatures determined for the doubly ionized nitrogen 
atom were not usable as they were consistently unrealistic and ir- 
regular. Temperatures varied from 10 electron volts to 237 electron 
volts, well above arcing potential. Evaluation of individual 
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temperatures showed that out of six temperatures determined for each 
machine run at constant magnetic field, four were negative temperatures. 
Investigation of the lines giving these data showed that each was an 
allowed transition, and that two of the lines are noted as among the 
more intense nitrogen lines. Improper identification of lines is not 
the problem, as in each case, there are no nitrogen lines or argon lines 
which could be misconstrued as the identified lines [4] . Their relative 
intensities are not consistent with those relative intensities shown 
in Table II, but the relative intensities of the singly ionized species 
are consistent. There is no readily acceptable explanation for these 
deviations because the relative intensities for all lines are determined 
in a single nan with a single base. As for the negative temperatures, 
it seems to follow that the doubly ionized nitrogen atcm is not in local 
thermal equilibrium, and that the velocity of the electrons is not Max- 
wellian. The electron density is probably such that radiative depop- 
ulation of lower energy states is prevalent, giving higher electoxn 
temperatures than the actual temperature. 

2. Relative Number Density, N++/N+ 

It was necessary following the results of the electron tempera- 
ture calculations to assume that the local electron temperature for the 
singly ionized nitrogen atom was the local electron temperature for the 
system for the purpose of obtaining relative number densities of the 
two species studied. The average local electron temperature for the 
singly ionized nitrogen atom for a particular magnetic field was entered 
as a constant in the relative number density program, and the correspond- 
ing computational data were entered and the program run. The data and 
the temperature were changed for each value of magnetic field studied 
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(See Figure 9). As the magnetic field increased to 4200 gauss, the 
relative number density decreased in accordance with the postulated 
increase in the density of the singly ionized atan. The density of 
the doubly ionized atom must also decrease somewhat, as the temperatures 
calculated show an approach towards the local equilibrium state, i.e., 
frcm 48 electron volts to 18 electron volts at 3600 gauss. This does 
not in any way assume that the condition of local thermal equilibrium 
is achieved by the doubly ionized nitrogen atcm. 

B. CONCLUSIONS 

This investigation has identified the variance in local electron 
temperature with magnetic field for the singly ionized nitrogen atcm 
and determined that the singly ionized nitrogen atcm is in local thermal 
equilibrium in the plasma column. The variance in temperature with mag- 
netic field of the doubly ionized nitrogen atcm was not determined to 
a satisfactory degree. The doubly ionized nitrogen atcm was determined 
not to be in local thermal equilibrium. 

It was shewn that the local electron temperature of the system was 
the local electron temperature of the singly ionized nitrogen atcm as 
evidenced by the comparison of local electron temperatures to relative 
number densities determined at the same temperature (Figure 10) . 

As the "theta pinch" device became operational late in the investig- 
ative period, and still requires further modification, it was not possible 
to include an analysis of shock wave perturbations on the nitrogen plasma 
column. It is felt that this can be accomplished in the future utilizing 
the equipment, as modified, and the procedures developed in this investi- 
gation . 
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ELECTRON TEMPERATURE COMPUTER PROGRAM 



THIS PROGRAM CALCULATES PLASMA ELECTRON TEMPERATURES FROM 
INTENSITY PATIOS ( GAI NEC RY PHOTOELECTRIC MEANS) OF SELECTED 
SPECTRAL LINES AMO THEIR CHARACTERISTIC CATA, ANC THEN, 

TAKES THESE TEMPERATURES AND CALCULATES THEIR MEAN, VARIANCE 
, AND STANDARD DEVIATION. 

DIME NS I ON GLOW(50),F(5O),EXCIPC(5O),RELINT(5C),WL(5O) 
lJLOW(EO) , JtJPPER (50) ,TEMP(500) ,RAT 10(500) 

DIMENSION ARMAX( 15) ,ARMIN( 15 ) 

PEAL GLOW, F.EXCIPC, RELINT, WL,0, TEMP, EVCCNV, RATIO, JLGW, 

IN, J UPPER, KRIS 
INTEGER S,Z 
VWL=500. F-l 0 
Vf =7.0 

THE ABOVE CARDS WILL PRINT CUT THE WAVELENGTH MAXIMUM ANC 
EXCITATION POTENTIAL MINIMUM DIFFERENCES USED IN THE CALCU- 
LATIONS. THFY SHOULD BE CHANGEC EVERY TIME THOSE DIFFERENCE 
S ARE CHANGED IN THE BODY OF THE PROGRAM. 

A P M AX ( 1 ) =6. 0 
ARMI N ( 11=1. 0/6. C 

THE ABOVE CARDS SET A MAX ANC A MIN ON THE THE RELATIVE 
INTENSITY PATIO OF ANY TWO (2) SPECTRAL LINES. 
EVCQNV=1.60207E-19 
B=l. 380<,2E-23 

1 FORMAT (15) 

2 FORMAT (4F10.0,E13.0,7X,F10.0) 

3 FORMAT ( '1 ' ,5X, 'NO. ' »4X , 'WAVELENGTH' , 7X» • EXCITATION 
1PCTENTIAL' »11X , 'J' ,17X, 'J' ,11X, 'ABSORPTION' » 4X , 

2'RE LATI VE',5X,'2J*1',/ 14 X,'(METERS)',15X, '(VOLTS)', 
311X, '(UPPER LEVEL )• ,5X, •( LOWER LEVEL )• ,6X, • STRENGTH ' , 
ASX, 'INTENSITY' ,4X , ' (LOW) '///) 

4 FCP M AT( /»5X,F4.C»2X,4PE13.3* 12X,0PF7. 3,17X,F3. 1, 1 5X, 
1F3. 1»11X,F8.5»6X»F7.2»5X,F4. 1 ) 

5 F T R M AT (• 1 • ,5X, 'NO. • ,5X, 'E-TEMPERATURES' ,5X, 

1 'WAVELENGTH! 1) • ,5X , • WAVELENGTH ( 2 )',5X,'INTENSITY RATIO 
2 • /13X, • ( DEGREES KELV I N ) • , 7X , ' ( M ET E PS ) • , 10X , • ( MET EPS ) • 
3///) 

6 FORMAT (5X,F4.0,6X,1PEU.4,7X,4PE13.3,5X,4FE13.3*8X* 

IIP Ell. A/) 

10 FORMAT <//////////25X, 'E-TEMPERATURE = • , 1 P E 1 1 . 4 , 2X , 

1 ' + /-', 2X,CPE11.4,2X, 'DEGREES KELVIN', 2X,' ************ ) 

11 FORMAT (2CX,4PE13.3/) 

12 FORMAT (• 1 • ,18X, • WAVELENGTHS CHECK'/////) 

14 FORMAT ( /////15X, 'THE TOTAL NUMBER OF ACCEPTABLE 

1CCNPUTATICNS (I. E., WITHIN SPECIFIED LIMITS) M A C E= ', 
2F5.0) 

200 FORMAT (• 1 • ,14X , 'MAXIMUM WAVELENGTH DIFFERENCE = ',3P£1 

IO. 1 ,2X, • METERS' ) 

201 FCRMAT( /////15X, • MINIMUM EXCITATION POTENTIAL 
1DIFFERFNCE = ' ,F4.1,2X, 'VOLTS'/////) 

202 FORMAT ( 15X, 'MAXIMLM ALLOWABLE INTENSITY RATIO = ', 
1F6. 1/////15X, • MINIMUM ALLOWABLE INTENSITY RATIO = ', 
2F5.2///) 

2 C 5 FORMAT ( // / //2 5 X , • E- TEMPERATURE = ' ,F7.4,2X, '♦/-• ,2X, 
1F7. 4, 2X, 'ELECTRON VOLTS' ,2X, •*****•) 

210 FORMAT ( //35X, • TEMPERATURE RANGE: • , 2X , 1 PE 1 1 . 4 , 2X , • TO • , 

12X,1PE1 1.4) 

211 FORMAT (//35X, 'TEMPERATURE RANGE: • , 2 X , F 7 . 4 , 2 X , • TO • , 

1 2X , F 7. 4 ) 

21 2 FORMAT ( //35X, ' PERCENT ERROR = ' ,2X,F7.2,2X,' PERCENT' ) 

2 C READ (5,1) J 
WRITE (6.12) 

IF (J) 100,100,21 
21 DC 80 1=1, J 

THE ABOVE DO LOOP CHECKS TO SEE IF ALL CATA HAS BEEN ENTEREC 
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READ (5.2) EXC I FO ( I) , JliPPER ( I ) , JLCW ( I) ,F ( I ) ,WL ( I ) , 
1RFLINK I ) 

WF ITE ( 6 • 1 1 ) WL ( I ) 
eC CONTINUE 

THE ABOVE INSTRUCTION WILL READ OATA CN A SERIES CF "J" 
SPECTRAL LINES. 

N = 1 . 0 

WRITE (6,3) 

DC 22 1=1, J 

GLCW< I ) =(2. 0*JLCW< I ) ) +1.0 

WRITE (6,4) N,WL( I ) ,EXC IPC( I ) , JUPPER( I ) , JLCW( I ) ,F( I ), 
IRE LI NT ( I ) , G L 0 W ( I > 

N=N + 1 . 0 

22 CCNTINUE 

DC 502 Z = 1 , 3 
AFMAX(Z+1)=ARMAX(Z)-1.0 
APMIN(7) = 1.0/ARVAX< Z + l) 

RATMAX=ARMAX ( Z + 1 ) 

R A TMI N=ARN I N ( Z ) 

N= 1 • 0 
S = 1 

WRITE (6,5) 

DC 26 1=1, J 
DC 66 MAP=I,J 

KPI S=EXCIPC(MAP)-EXCIFC( I ) 

CHECK=ABS(KRIS) 

IF (CHECK. LT. VE )GC TC 70 

THE ABCVE INSTRUCTION ChECKS TC SEE IF THE DIFFERENCE IN EX- 
CITATION POTENTIAL IS GREATER THAN 7.0 ELECTRON VOLTS. THE 
DIFFERENCE MUST BE SEVERAL VCLTS FCR VALID RESULTS. 

WAVF= ABS ( WL ( I ) -WL ( HAP ) ) 

IF (W'AVE.GT. VWOGC TC 70 

THE ABOVE INSTRUCTION CHECKS TO SEE IF THE DIFFERENCE IN 
WAVELENGTHS IS LESS THAN 500 ANGSTRC^S. MUST BE FCR VALID 
RESLLTS. 

IF (KRIS) 2 4,100,23 

23 PATIO(S)=RFLINT(I) /RELINT (MAP) 

VANG=PATI 0 ( S ) 

IF (VANG.GT.RATMAX)GC TC 70 
IF ( VANG. L T. RATH IN ) GC TC 70 

TEMP(S)=(((EXCIPO(MAP)-EXCIPC(I) >*<EVCONV) )/B)/ALCG 
1( (PATIO(S) ) * ( G LC W ( MAP)/GLOW( I ) )*(F(MAP)/F( I) ) * ( WL ( I ) ** 
23.0/WL( MAP) * *3 . C ) ) 

GO TO 25 

24 RATIC(S)=RELINT(MAP)/RELINT(I) 

VUNG = R AT I C ( S ) 

IF ( VUNG.GT. RATMAXJGC TC 70 
IF ( VUNG.LT. RATMI N ) GC TC 70 

TEMP(S) = ( ( (EXCIPO( I J-EXCIPC(MAP) )*(EVCCNV) )/B J/ALCG 
1( (PATIO(S) ) * ( GLCW ( I )/GLCW(MAP) ) * ( F ( I )/F(MAP) )*(WL (MAP ) 
2**3 . C /WL ( I ) **3 . C ) ) 

GC TO 25 
7 C TEMP(S)=0.0 
RATI C( S ) =0. 0 

THIS PRL'GP AM IS SET LF TC RECCRC A VALUE CF 0.0 FCR TEMPERA- 
TURE AND RATIO IF THE EXCITATICN POTENTIAL DIFFERENCE OR 
WAVELENGTH DIFFERENCE DC NCT MEET ABOVE CRITERIA. IT WILL 
ALSO ADJUST THE VALUE CF N FOR THE COMPUTATION CF THE MEAN, 
ETC. . 

25 WRITE (6,6) N , TE M F ( S ) , WL ( I ) , WL ( M A P ) , R AT 10 ( S ) 

N = N+ 1 • 0 

S = S+1 

66 CONTINUE 

26 CONTINUE 
SL M =0. 0 
ASLM=C. 0 
N = N- 1 

S = S-1 

DC 27 1=1, S 

TEMP ( I )=AeS( TE M F ( I ) ) 

IF ( TE MP ( I ) ,E0.C.C)N = N-1.0 

27 SLM=SUM+TEMP( I ) 
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IF (N.LE.l.O)GC TC 5C2 
AVF=SUM/N 

IF (N.EQ.l.O)GC TC 5C2 
DC 28 1=1, S 

IF ( TEMP( I) .EQ. 0.0)GC TC 123 
PIFF=A8S( AVE-TEMP( I ) ) 
SC=DTFF**2 
Gf TC 128 

123 DIFF=0.0 
SC=0.0 

126 A SLM = A SUM + S G 

S I GM= SORT I A SUM/ (N- 1.0) ) 
SIGP=SIGM/SGRT(N) 

WRITE (6,200) VWL 

WRITE (6,201) V E 

WFITfc (6,202) RATMAX ,RATMIN 

WRITE (6,14) N 

WRITE (6,10) AVE,SIGR 

P A VE=A VE + S I GB 

0/ VE=AVE-S!GB 

WFITF (6,210) PAVE, CAVE 

PERR=SIGP*100.0/AVE 

WRITE (6,212) FERR 

EAVE=AVE*P/EVCCNV 

E$IGP=SICB*B/EVCCNV 

WRITE (6,205) EAVE,ESIGB 

PEAVE=EAVE+SIGP 

CEAVE=EAVE-SI Gfi 

WFITE (6,211) REAVE, CEAVE 

502 CONTINUE 
GC TO 20 

ICC STOP 
END 
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RELATIVE NUMBER DENSITY COMPUTER PROGRAM 



THIS PROGRAM CALCULATES THE RELATIVE NUMBER CENSITIES OF TWO 
SPECIES IN A PLASMA FRC M THEIR SPECTRAL LINE DATA. THE TWO 
SPECIES APE CALLED NL AND N2.WITH THE N2 SPECIES HAVING THE 
HIGHER EXCITATION POTENTIAL (IN GENERAL) AND BEING ENTERED 
AS THE SECOND OF TWO CATA SETS. 

DIMENSION GLOW (4, 50) f F(4 ,50) , EXCIP0(4,50) f REL INT(4,50) 
1 WL ( 4 , 50 ) , JL0W(4,5 0) , JUPP ER ( 4 , 50 ) , R AT I 0 ( 50 , 50 ) , 

2NUMDEN ( 50,50) 

REAL GLOW.F .EXCIPC ,RELINT,B,WL, TEMP, EVCON V , R AT 1 0 , JLOW, 
1JUPPER , KRIS, NONE, NTWC.NUMDEN.N, MATE 



EVCONV=1.30207E-1'9 
B=l. 38042E-23 
VWL=500 . E-l 0 
VE =7. 0 

THE ABOVE 2 CARDS GIVE THE MAXIMUM WAVELENGTH AND MINIMUM 
EXCITATION POTENTIAL DIFFERENCES USED IN THE COMPUTATIONS. 

R ATMA X=4. 0 
R A TMI N=0. 2 5 

THE TWO CARDS ABOVE SET A MAX AND A MIN CN THE RELATIVE 
INTENSITY RATIO OF ANY TWC SPECTRAL LINES. 

TEMP=44700. 0 

TEMPERATURE IS CONSIDERED A CONSTANT FOR THESE CALCULATIONS. 
THE VALUE OF TEMPERATURE USED SHOULD BE OBTAINED FRCM THE 
ELECTRON TEMPERATURE CCMFUTER PROGRAM PGR THE SPECTRAL CATA 
USED HFREIN. 

1 FORMAT! * 1 * ,9X, ’ELECTRON TEMPERATURE = ',F9.0,2X, 

1 'DEGREES KELVIN*) 

2 FORMAT (15) 



FORMAT ( • 1 • ,20X, • CATA CHECKS*/////) 
FORMAT (4F10.0 , E 1 3 . 0 ,7X , F 10 . 0 ) 
FORMAT ( 2 OX ,4PE 1 3.3/) 

FORMAT ( *1 * ,3X , *N0. • ,5X, • WAVELENGTH 



_ ,5X, 'WAVELENGTH* ,7X. • EXCITATION 

1PCTENTIAL*,11X,*J*,17X,*J*,11X, • ABSORPT ION • ,4X. 

2* RELATIVE • , 5X , *2 J* I • / 13X . • ( METERS ) • , 1 5X, • ( VCLTS ) • , 1 IX , 
3 • ( UPP E R LEVEL) *,5X, ‘(LOWER LEVEL )*, 6X ,* STR ENGTH *, 5X, 

4* INTENSITY' ,4X, ' ( LCW) •///) 

7 FORMAT ( 1X,F4.0,F4.C,IX,4PEI3.3,I2X,0PF7.3,I6X,F4.1, 
1I4X.F4. 1,I0X,F8.5,6X,F7.2,5X,F4.1/ ) 

8 FORMAT ( *1 • ,4X , *N2. • , 3X, • WAVEL ENGTH ( 2 ) • ,5X,*N1. • ,3X, 

1 • WAVELENGTH! 1) • , I OX , • I NT ENS I T Y RAT 1 0 • , 1 1 X , • NUMBER 
2DENSITY RATI OS '/l 3X ,* (METERS )*, I6X ,*( METERS )', 15X, 
3'(I(2)/I(1))',19X,'N(2)/N(1)'///) 

<3 FORMAT ( /4X,F4.0,3X,4PEI3.3,4X,0PF4.0,3X,4PEI3.3, I2X, 
LIFE 10. 3,20X,2PE10.2 ) 

10 FORMAT ( • 1 • ,14X, • SUM = • , 4P E 12 . 2 , / / 1 5X , • A VER AGE 

INUMBER DENSITY RATIC = • , 2 PE 10 . 2/ / / / / 1 5 X , • A 8 S ( NB AR-N ) 

2* ,15X, 'ABS( (NBAP-N)SO) ' ) 

11 FORMAT ( /15X.2PE10.2 , 18 X , 2PE 10 . 2 ) 

12 FORMAT ( • 1 • ,14X, • THE TCTAL NUMBER OF ACCEPTABLE 

1C CMPUT AT I CNS (I. E., WITHIN SPECIFIED LIMITS) MACE = *, 

2F5.0) 



13 F0RMAT(///15X, 'STANDARD 
1, 'STANDARD DEVIATION CF 



CEVIATICN = 
THE M EAN * 



• ,12X,EIC.4//15X 

' , E 1 0 • 4 ) 



14 FORMAT ( ///////////20X,'THE NUMBER DENSITY RATIO OF 
1SPECIES 2 TO SPECIES 1 = • ,2PE 10.2 , 2X, •«■/-', IX, 0PE8. 2 

200 FORMAT ( / / / / /I OX , • M A X I MUM WAVELENGTH DIFFERENCE = ', 

13FE10. 1, 2X, 'METERS' ) 

201 FORMAT (/////10X, 'MINIMUM EXCITATION POTENTIAL 

1DIFFFFERENCE = ', F4. 1 , 2X ,' VCLTS •///// ) 

202 FORMAT ( 10X , • MAX I MUM ALLOWABLE INTENSITY RATIC = ', 

1F6.1/////10X, 'MINIMUM ALLOWABLE INTENSITY RATIC = ', 
2F5.2///J 

SPECIES 2 = N++ , AND SPECIES 1 = N+ FOR THIS INVESTIGATION. 
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5CC RE AD ( 5.2 ) INOIA 

IF (INDIA) 501,501,502 
5C2 WRITE (6,3) 

DP 30 1=1, INDIA 

READ (5, A) EXCI F0( 1 ,1 ) , JUPPER( 1 , I ) , JLCW( 1 , I ) ,F( 1, I ) , 
1WL(1,I) .RElINT(l.I) 

THE ABOVE CARD WILL R E AC THE M DATA SET CONSISTING OF A 
TOTAL OF "INDIA” NUMBER OF CARCS. 

WRITE (6,5) WL( 1,1 ) 

THIS PRINT-OUT CHECKS TC SEE THAT ALL CATA IS ENTEREC. 

30 CONTINUE 

READ (5,2) JULIET 
WRITE (6.3) 

DC 31 1=1, JULIET 

READ (5, A) EXCIP0(2»I ) »JUPPER(2» I ),JL0W(2» I),F(2, I) , 
1WL (2,1 ) ,PELINT( 2, I ) 

THE ABOVE CARD WILL REAC THE N2 DATA SET CONSISTING OF A 
TCTAL OF "JULIET" NUMBER OF CARDS. 

WRITE (6,5) WL ( 2 , I ) 

THIS PRINT-OUT CHECKS TC SEE THAT ALL DATA IS ENTEREC. 

31 CONTINUE 
NO N'F = 1 • 0 
N T W0= 1.0 
WRITE (6,6) 

K I LP= I ND I A 
DC 33 K=1 ,2 

DP 32 I = 1 , K I LC 

GLPW(K.I)=(2.0*JLCW(K,I))+1.0 

WRITE (6,7) NONE,NTWC,WL(K,I ) , EXC I PC ( K , I ) , JUP P ER ( K , I ) , 
1JL0W(K,I),F(K,I),RELINT(K,I),GLCW(K,I) 

FCR THE RESULTING PRINT-OUT, THE FIRST NUMBER WILL INDICATE 
THE DATA SET AND THE SECOND WHICH ONE OF THAT DATA SET. 

NT WO=NT W0+ 1.0 

32 CONTINUE 
NCNE=NONE+ 1.0 
N TWO= 1.0 

K I L0 = J UL I E T 

33 CONTINUE 
NONE = 1 . 0 
NTW0=1 . 0 
WRITE (6,8) 

DC 39 K=l, JULIET 
DC 38 1=1 .INDIA 
KPIS=EXCIPP< 2,K)-EXCI PC (1,1 ) 

CHECK=ABS(KRIS) 

IF ( CHECK. LT.VE)GC TC 3A 

THIS INSURES THAT THE DIFFERENCE IN EXCITATION POTENTIALS IS 
GREATER THAN OP EQUAL TC 7.0 VOLTS WHICH IS NECESSARY FCR 
VALID RESULTS. 

MATE=AB$(WL(2,K)-WL( 1 , 1 ) ) 

IF ( MATE. GT. VWL ) GC TC 34 

THIS INSURES THAT THE DIFFERENCE IN WAVELENGTHS IS EQUAL TO 
OR LESS THAN 500 ANGSTROMS. NECESSARY FCR VALIC RESULTS. 
RATIC(K f I )=REL I NT (2 ,K ) /RELINT( 1 , I ) 

VCNG=RATI 0 ( K , I ) 

IF ( VCNG.GT. RATMAX ) GC TC 34 
IF (VONG.LT.RATMIN)GC TC 34 

NUMDEN ( K , I ) = ( RATI C(K, I )*(GLOW( 1,1 )/GL0W(2»K))*(F( 1. I)/ 
1F(2,K) )*C ( WL(2 ,K) /WL( 1 , I ) )**3.0 )*(EXP( (KRIS*EVCONV)/ 

2 ( B*TEMP ) ) ) ) 

GC TP 35 

34 NUMDE N ( K • I )=0. 0 
RATI 0( K ,f ) =0.0 

IF THE EXCITATION POTENTIAL DIFFERENCE CR WAVELENGTH DIFFER- 
ENCE DO NOT MEET ABOVE CRITERIA, THE VALUES OF THE NUMBER 
DENSITY RATIO ANC THE INTENSITY RATIO WILL eE SET TC ZERO. 

35 WFITE (6,9)NTWC,WL(2,K),N0NE,WL(1,I),RATIC(K,I), 

1 NUMPE N ( K , I ) 

NCNF = NONE+ 1.0 
XP A Y=NCNE 
38 CONTINUE 

NT WO=NTWP+ 1.0 
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ZFBRA=NTWC 
NONF = 1.0 
3 5 CCNTINUE 

N= ( XRA Y- 1 . 0 ) * ( ZEBRA- 1.0) 

SUM=0.0 
ASUM=0. 0 

DC 46 K = 1 , JUL I E T 

DO 45 1=1. INDIA 

NL'MDEN(K,I ) =ABS(NUMDEN(K, I ) ) 

IF (NUMDEN(K.I ).EC.0.0)N=N-1.0 
THE ABOVE CHECK CORRECTS THE COUNT TO. ACCEPTABLE CCMPUTA- 
TIONS ONLY, IN CPDER TO KEEP. THE STATISTICAL DISTRIBUTION 
C CR R EC T • 

4 5 SUM=SUM+NLNDEN(K,I ) 

46 CCNTINUE 
AVE=SUM/N 

WRITE (6,10) SUP , A VE 
DO 48 K=l, JULIET 
DC 47, 1=1, INDIA 

IF (NUMD£N(K,I ).EC.O.O)GC TO 123 
DIFF=ABS( AVE-NUPDEN(K, I ) ) 

S0=DIFF*2 
GC TO 124 

123 DIFF=0.0 
S0 = 0.0 

124 WRITE (6,11) D I FF , SC 

47 ASUM=ASUM+SQ 

48 CONTINUE 

SIGM=SQPT( ASUM/( N-l .0) ) 

SIGB=SIGM/SQRT(M 

WRITE (6,12) N 

WRITE (6,13) SIGP,SIGB 

WRITE (6,1) TEPF 

WRITE (6,200) VWL 

WRITE (6,201) VE 

WRITE (6,202) RATPAX , RATPI N 

WRITE (6,1M AVE , SIGR 

THE PRINT-OUT ABOVE GIVES THE AVERAGE NUMBER DENSITY (N2/N1) 
+ /- THE STANDARD DEVIATION CF THE MEAN. 

5 C 1 STOP 
END 



SAMPLE DATA CARDS 



(NCT TO 



SCALE ) : 



E 

g 


J(U) 


J ( L ) 


F 


WAVELENGTH 


I (REL 


28.360 


2. 0 


1.0 


.118 


+4124.100E-1C 


21.49 


28.361 


2. 0 


2.0 


.118 


+4133.654E-10 


3 5.07 


28.362 


2. 0 


3 . C 


.118 


+4145.759E-1C 


41.52 


21.070 


2.0 


1.0 


.133 


+4601 . 490E-1 0 


96.65 


21.062 


1.0 


0.0 


.320 


+4607. 167E-1C 


89.90 


21.061 


1.0 


1.0 


.0798 


+4613.887E-1C 


60.23 


21.060 


0. 0 


1. C 


. 106 


+4621.405E-1C 


89.44 


21. 06* 


1.0 


2.0 


.0796 


+4643. 106E-10 


107.50 


g 

3C. 331 


1.5 


0. 5 


.484 


+4097.309E-1C 


13.74 


3C. 330 


0. 5 


0.5 


.242 


+ 4103.3 70E- 1 0 


10.46 


35.640 


2.5 


1.5 


.398 


+4200. OOOE- 1 0 


157.05 


38.250 


3. 5 


2. 5 


.286 


+4514 . 900E- 1 0 


12.68 


38.230 


0. 5 


0.5 


. 178 


+4518.200E-1C 


52.06 


41.510 

-ICO 


1.5 


1.5 


. 304 


+4535. 100E-1C 


22.04 
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